Analysis of volume distribution of power loss in ferrite cores 

M. LoBue,^ V. Loyau/ and F. Mazaleyrat^ 

SATIE, ENS de Cachan, CNRS, UniverSud, 61 av du President Wilson, 
F-94230 Cachan, France 

(Dated: 18 January 2013) 

We present a technique to estimate the inhomogeneities of magnetic loss across the 
section of ferrite cores under AC excitation. The technique is based on two distinct 
calorimetric methods that we presented elsewhere. Both the method are based on 
the measurement of the rate of increase of the sample temperature under adiabatic 
condition. The temperature ramp is recorded either measuring the sample bulk 
resistivity, or using a platinum probe pasted on the sample surface. As an example 
we apply the procedure to an industrial sample of Mn-Zn ferrite under controlled 
sinusoidal excitation with a peak induction of 50 mT in the range between lOOkHz and 
2MHz. The results are discussed by comparison with simulations of the dissipation 
field profile through the sample, calculated using a FEM code. 
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I. INTRODUCTION 



Soft sintered Mn-Zn ferrites are widely used in high frequency apphcations based on 
switching circuits. Their relatively high resistivity and permeability make them optimal for 
the 30kHz- IMHz bandwidth. Miniaturization and the increase of working frequencies make 
thermal and loss managements the keys for future design of electronic switching devices. It 
is thus necessary to establish reliable methods to measure, understand and predict the losses 
in soft ferrites at medium and high frequencies. 

Losses are commonly measured using the flux-metric method^"-. However, when excita- 
tion frequencies are above 10 kHz, measurements can be dramatically affected by spurious 
phase contributions due to stray capacitances and inductances of the circuits and probes. A 
possible alternative is the calorimetric method; in adiabatic conditions, the power dissipated 
in the magnetic material is directly related to the heating rate, reducing the measurement 
to that of a temperature ramp. We discussed in detail this technique in a recent paper- 
where an experimental set-up using a platinum probe as temperature sensor was presented. 
A further improvement of the technique has been proposed^ by measuring the temperature 
from the bulk DC resistivity of the ferrite core^. Whereas the temperature deduced from 
bulk resistivity represents a measurement averaged over the whole sample volume (as in the 
case of flux- metric measurements), the platinum probe measurement gives information on 
dissipation taking place in a thin strip of the material just near its surface. Consequently 
the comparison between the two measurements can be used to estimate the loss inhomo- 
geneity, and thus, the relevance of skin effect for a given sample geometry and excitation 
frequency. It will be shown that, global low level electromagnetic parameters, obtained by 
impedance spectroscopy, can be used as local input parameters for the FEM computation of 
magnetic loss in linear regime. In this way we shall calculate the volume or surface averaged 
loss values, corresponding to the two types of calorimetric measurements and compare them 
with experimental results. 



II. EXPERIMENTAL AND COMPUTATION DETAILS 



Power losses in magnetic materials under AC excitation are associated with all the ir- 
reversible phenomena taking place during the core magnetisation (i.e. Barkhausen jumps, 
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eddy currents, domain wall resonance, spin dumping, magnetoacoustic emission, etc.). Un- 
der adiabatic conditions^, all the dissipated power is converted into heat and induces a tem- 
perature increase of the sample. Thus the temperature field through the core T{x,y, z,t) 
can be described by means of the linearised heat equation giving: 

dT{x,y,z,t) 1 

where ps{x,y,z) is the local heat source density (in Wm~^), Cp is the specific heat (in 
JK~^kg"^) of the sample and pm its density (in kgm~^). Now, if Vq is the volume of the 
sample, the magnetic loss per cycle measured using the fiux-metric method, at a given 
excitation frequency /, corresponds to the average loss over the sample volume, namely: 

iW) = ^ I ^^(fl^rfV. (2) 

Whereas a measure of the average temperature of the whole sample will result in a correct 
extrapolation of (W), any local measurement will be related with dissipations taking place 
just in a small volume near the position of the temperature probe. In the latter case (yV) 
can be deduced properly only when volume losses are uniformly distributed through the 
sample. 

Two different techniques where used to measure the temperature ramp: the first using 
a platinum probe pasted (PtSOO) on the sample^ and therefore sensible to the temperature 
change in a narrow strip just near the surface, and a second technique based on the mea- 
surement of previously calibrated bulk resistivity of the ferrite core- which returns a global 
average of the temperature ramp. The latter leads iyV) notwithstanding the inhomogeneities 
of the loss field due to the geometry of the sample and to the skin effect. In both cases, 
the sample was a square cross section bar (I rod: 0.64 x 0.64 x 2.54 cm^) made of 3E27 
Ferroxcube Mn-Zn ferrite. Complex impedance spectra of the permeability /i = fi' + ifi" and 
of the electrical conductivity a = a' + ia" = a' + itu (where e is the electrical permittivity 
and UJ = 27r/ is the angular frequency) where performed with a HP 4195 impedance analyser 
between 100 kHz and 10 MHz. The losses where measured in a vacuum chamber between 
100 kHz and 2 MHz as described in ref.->^ from 10 second long temperature recordings. 

In order to validate the interpretation of the measurements performed with the two 
techniques, an approach based on the solution of Maxwell equations was used. For this 
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FIG. 1. Total loss field wtot{x,y) calculated from Eq. (g]) and ([5]) for / = lOOkHz [(a)] and f=2MHz 
|(b)[ The dashed rectangle indicates the region used to calculate the surface loss as measured by the 
platinum probe. The thickness of the rectangle is defined as the product of the heat propagation 
speed in the material and the measure duration (10 seconds). Grey scales indicate the local loss 
in mJkg~^ where the loss varies from 0.04mJkg~^ (dark regions) to 0.05mJkg~^ (white regions) 
(a) , and from 0.45mJkg^^ to 0.85mJkg^^ in (b) . The ratio between the maximum and the 
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minimum values is 1.25 at / = lOOkHz [(a)] and 2 at f=2MHz [(b)] 



purpose, the code FEMM^ was used taking measurements of complex magnetic permeability 
and conductivity as intrinsic material parameters. All the following analysis is based on 
three strong assumptions: {i) the peak induction values are sufficiently low to consider the 
magnetic behaviour as linear. Thus the concept of real and imaginary magnetic permeability 
and their relationship with energy losses is valid. This assumption is confirmed by the fact 
that all the loss measurements presented in this paper scale roughly as the square of the peak 
induction. This means that the term /i"/ | P remains constant in the investigated induction 
range (from 3mT to 50mT); (m) all the measurements of complex magnetic permeability are 
performed on a toroidal ferrite sample with section small enough to make eddy current loss 
and skin effect negligible; {in) all the measurements of complex conductivity are performed 
on bar shaped samples which section is sufficiently reduced to neglect skin effect. 

From (m) the measured /x is a parameter related to the local phase lag between H and B 
due to all dissipative phenomena but eddy current losses. On the other hand, the complex 
electrical conductivity a results from the mixed electrical properties of the grain and of the 
grain boundaries^. Its real part cr', is related to ohmic and dielectric loss effect whereas the 
imaginary part a", is due to the effect of mixed permittivity. 
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FIG. 2. Comparison between the bulk measurements of the average total loss and the calculated 
ones. Peak inductions are: 3mT, lOmT, 30mT, 40mT and 50mT 
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FIG. 3. Comparison between bulk and surface loss measurements and the same for computations 
(over the whole volume and the one delimited by the dashed line in FigHJ The peak induction is 
of 50mT 

Given the geometry of the problem, with the magnetic field always perpendicular (along 
z axis) to the core section, the calculation of the loss field reduces to a two dimensional 
problem in the xy plane. Now, under sinusoidal excitation when Hz{x,y,t) = HQ{x,y)e^'^^, 
substituting the measured a and /i in the diffusion equation for the magnetic field we have. 



- V^i^^(x, y, t) = iuaHo{x, y). (3) 
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From this equation and from the fact that j = V x H we can calculate the eddy-current 
loss field across the core section, 



Pecix,y) = -p\j{x,y)\^ (4) 

where p' = Reil/a). This expression includes all the eddy current related losses. The 
effect of ohmic and dielectric dissipation in both, the grains and their boundaries, are de- 
scribed by the real part of the resistivity p' . On the other hand, the complex part of the 
permeability /i", measured on the very small sample, will account for all the magnetisation 
process related dissipation phenomena ( excess losses, spin damping, etc.). So we can write 
the magnetic loss field as: 

1 / I p" 

Pmagix,y) = -uj p" \U,{x , y)\^ = -uj-^\B,{x , y)\^ (5) 

The total power loss will be the sum of the two terms of Eq. and ([5]), and the loss per 
cycle and unit mass will be Wtot{x,y) = Ptot{x,y) /{f pm)- By averaging over the core section 
and multiplying by the length, we obtain (W)^ = i /"^^2 /-f/2 "^(^i y)dxdy corresponding to 
the bulk (resistance) measurement, whereas {W)^ = £ J^^'^ w{x,y)dxdy yields the value 
corresponding to surface (PT 500) measurement. The thickness d is defined as the product of 
the heat diffusion velocity in the material and the duration of the measurement (10 seconds). 



III. RESULTS AND DISCUSSION 

In order to determine the magnetic field Hz{x, y) across the section area we must specify 
the induction conditions under which losses have been both calculated and measured. Here 
we consider the case where a sinusoidal magnetic induction of average value (Bq) is forced 
across the core section area S. This implies that the magnetic field verifies: 

{Bo) = l f pHo{x,y)dS (6) 
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S 

The magnetic field Hz{x, y) has been determined solving Eq.([3]) coupled with Eq.([n]) with 
a two dimensional finite element code in a frequency range between lOOkHz and 2MHz on 
the sample geometry. Calculation results are shown in Fig. [1] where the loss field Wtot{x,y) 



is plotted for / = lOOkHz in Fig 1(a) and / = 2MHz in Fig, 1(b) 



In Figj2] comparison between the bulk measurements of loss and the calculated average 
total loss (W)^ is shown for peak inductions ranging from 3mT to 50mT. The agreement be- 
tween measurements and calculations is particularly good for frequencies above 300kHz. The 
over estimation of computed losses may be due to the compensation error in the measure- 
ment of fi" performed with the impedance analyser. Indeed, at low frequency the resistive 
part of the impedance is nearly zero, so even a small absolute error can yield a strong system- 
atic error. Fig l3] shows the comparison between the surface measured loss and corresponding 
computation (W)^. Calculations and measurements are in good agreement. With the sam- 
ple geometry studied here, as the average magnetic path is not much different from the 
extremal paths in the sample, all the loss inhomogeneities can be ascribed to the skin effect. 
So, either measurements and finite elements calculations, show that skin effect is negligible in 
a frequency regime below 600-700kHz. Above these frequencies skin effect becomes relevant, 
and a large part of the dissipation phenomena is concentrated near the surface of the sample. 
In conclusion, it has been shown that the calorimetric loss measurement method, exploiting 
bulk and surface temperature recordings, is able to give indirect information on the eddy 
current losses in Mn-Zn ferrites. Using FEM computation -for which material parameters 
input where obtained by impedance spectroscopy- it is possible to compute the loss fields 
and to obtain very good agreement with measurements without introducing any fitting pa- 
rameter. Hence, this model could be integrated in magnetic components design software in 
order take into account the magnetic losses including eddy currents for the dimensioning. 
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